ESTIMATES ON THE SPECTRUM OF FRACTALS 
ARISING FROM AFFINE ITERATIONS 



Palle E.T. Jorgensen and Steen Pedersen 

Abstract. In the first section we review recent results on the harmonic analysis 
of fractals generated by iterated function systems with emphasis on spectral duality. 
Classical harmonic analysis is typically based on groups whereas the fractals are most 
often not groups. We show that nonetheless those fractals that come from iteration of 
affine mappings in R d have a spectral duality which is instead based on approximation 
and a certain dual affine system on the Fourier transform side. The present work is 
based on iteration of frame estimates (which have been studied earlier for regions in 
K d ). Our emphasis is on new results regarding the interplay between the limit-fractal 
on the one hand, and on the other the corresponding regions in M d which generate 
iterated function systems of contractive affine mappings. As an application of our 
frame results, we obtain a classification of a certain type of spectral pairs. 



1. Introduction 

Let f2 be a measurable subset of M. d {d — 1,2,...) of finite positive measure, 
i.e., for fixed d, denotes the (i-dimensional Lebesgue measure and we assume, 
< m d (il) < oo. Then the Hilbert space C 2 (il) is formed relative to the usual 
inner product 

(1-1) (f,g)n := / W)g{x)dx 

Jn 

where dx :— dxi . . . dxd, and ||/||q : = (/, f)n- If A = (Ai, . . . , A^) G M. d , then we 
shall need the exponential functions 

(1.2) e x {x) := e i2nX ' x = e a7r(A 1 x 1 +-+A <i x <J )_ 

These functions are defined on all of M. d of course, but, when f2 is given as above, 
they restrict and define elements in £ 2 (Q), viz., xne\ where xn denotes the indi- 
cator function of f2. It will be convenient to denote these restricted functions also 
by eA- If A is a subset of R d such that the corresponding exponential functions 
{eA : A € A} form an orthogonal basis for the Hilbert space £ 2 (£1), then we say 
that (57, A) is a spectral pair. Such pairs have been studied extensively recently, see 
e.g., [Fu], [Jol], [JoPe3-4], and [LaWa]. Defining the Fourier transform 

(1.3) /(A) := / V@jf(x) dx = (e x J) n 
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we then get 

(1.4) £|/(A)| 2 = m d (ft) / \f(x)\ 2 dx = m d (0)\\f\\ 2 a . 

Of course, the set f2 is implicit in the notation on the left-hand side of (1.4), as O 
enters into the definition of the transform, / /, see (1.3). 

It turns out that the structure of spectral pairs is surprisingly rigid, and it is also 
related to tiling questions, see [Fu], [JoPe3] and [LaWa]. The generalized spectral 
pairs (see details below) turn out to form a much richer class. They arise from a 
relaxation of (1.4), transforming it instead into a frame estimate, see also [Benl-2, 
BeTo, Ga, DuSc, Beu]. We say that (f2, A) is a generalized spectral pair (g.s.p.) 
in R d if f2 is as specified, and if there are constants, k, K, positive and finite, but 
depending on both and A, such that 

(i-5) wii£<£i/(a)i 2 < wii^ 

AgA 

It can be checked that if A is a set satisfying (1.5) for as specified, then A must 
be discrete. 

Lemma 1.1. Let (fi,A) be a generalized spectral pair, and let 
K = A := {t G R d : t ■ A G Z, VA G A}. 



Then K is discrete. 

Proof. The set {e^ : A G A} is total in £ 2 (£1). Hence by [Fu], A cannot be contained 
in any hyperplane in M. d ; in particular, there exists a basis {vi, . . . , v^} for WL d with 
each vj G A. Let w\, . . . , wa be a dual basis, viz., wj ■ Vi = Sji, for i, j = 1, . . . , d. 
Then we get the containment, 

K c {t E R d : t ■ v 3 ■ e Z, j = 1, . . . , d} 
= ^YltjWj : tj e z| , 

with the latter set a lattice in M. d ; hence K is discrete. □ 

Note, we only need {e\ : A G A} to be total in £ 2 (S!) for A to be discrete — not 
any of the other ingredients in the definition of a generalized spectral pair are in 
fact needed for this part. 

It is also know that, for the spectral pairs, if G A, then the group 

(1.6) A°:={sGR d :s-AGZ, VAgA} 

is of rank d; in other words, A is a lattice in M. d . If A is also a lattice, we have 
just the familiar situation of multivariable Fourier transforms; but there are many 
extremely interesting cases where (fi, A) is a spectral pair, but A is not a lattice. We 
take up this question in Sections 2.3-2.4 below which provides in fact a classification 
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of some of the possibilities. Although the results there are about spectral pairs, the 
arguments use g.s.p.'s and iteration. 

If (ft, A) is a generalized spectral pair, then there is a well defined operator 
F := F (n;A) from C 2 (Q) into £ 2 (A), given by 



(1.7) 



(F/)(A) := /(A), A e A 



and it is called the associated frame operator. We say that the system is exact if 
F(n.A) maps onto £ 2 (A), and that it is tight if k = K in (1.5). It follows from 
Hilbert space theory that the system is exact if and only if there are no non-zero 
sequence solutions (£a)asA m ^ 2 (A) such that the corresponding function f(x) — 
YlxeA^ e >-( x ) 1S zer0 as an element in £ 2 (f2). 

Let R be a real d by d matrix which is expansive; i.e., we assume that the 
eigenvalues Xi of R all satisfy \\\ > 1. If R = (r^), the corresponding transpose 
(rji) will be denoted R* . Two systems of affine maps will be considered as follows: 
Let B and L be given finite subsets in R d , and set 



(1.8) 
(1.9) 

(1.10) 



a b (x) 

n(s) 



R- 1 (x + b)=R- 1 x + R- 1 b, 
R*(s + £) = R*s + R*£, 

U °*( fi )> 

b£B 



and 



(1.11) r(A) := |J n(A). 

leL 

The affine mappings (1.8)— (1.9) are defined for x, s 6 M. d , and they are indexed by 
the respective (finite) sets B and L. In (1.8)— (1.11), we are applying these maps to 
subsets of R d , e.g., the components of some pair (fi, A) as above, or to more general 
sets. We shall be primarily interested in the case when the unions in (1.10)— (1.11) 
are non-overlapping, i.e., when 



(1.12) a b (Q) nov(Q) = for all b ^ b' in B 
and 

(1.13) n(A) riTp(A) = for all £ ^ £' in L. 



As far as (1.12) is concerned, we will relax it and allow overlap on sets of d- 
dimcnsional Lebesgue measure zero, i.e. 

(1.14) m d (a b (Q) nov(Q)) = for all b ^ b' in B. 

But it turns out that the harmonic analysis of systems subjected to (1.14) may be 
reduced to that of an associated system which in fact satisfies the stronger property 
(1.12), so we will restrict here to the latter case. When (1.14) holds we say that 
the triple (R,Q,B) is an affine system with no overlap, and if (1.13) holds that 



4 



PALLE E.T. JORGENSEN AND STEEN PEDERSEN 



(R* , A, L) is a dwaZ affine system with non-overlap. Starting with (1.10)— (1.11), wc 
define recursively, flo ■= 0, A := A, 

(1.15) n n+1 := a(n n ), 

and 

(1.16) A„ +1 := r(A„). 

For each n, let /x„ be the normalized measure coming from restricting Lebesgue 

measure to fi„, i.e., 

(1-17) /X„ : = Xf2„ X TOdl^n)" 1 X 77J d . 

The interpretation of (1.17) is 

/x„(A) = md(f2 n ) _1 md(A fl fi„) for all Borel sets A, 
or equivalently, 

/ / d , n = m d iSl n r J fXn n dm d for all Borel functions /. 
It can be checked that then 

(!-!8) Mn+i = 757 Yl ° o-j 

where -B denotes the number of elements in the (finite) set -B, and ix„ o er^ 1 is the 
measure given by, 

fj n oa^ 1 (E) := ^(o^ 1 ^)), 

where 

o- 6 " 1 (£;) = {x:ai(a;) ££}, 

or, equivalently, 

(1-19) /" / (d/Xn ° ^fc" 1 ) = /(/ O <J b ) dn n 

for Borel subsets E, and Borel functions /. Using then a standard theorem (see 
[Hut] and [Fal]) there is a unique limit measure /x (i.e., /x n — > /x in the Hutchinson- 
Hausdorff metric), and 

(1-20) H= T^E^ ^" 1 - 

This is typically a fractal measure, and the object of the present paper is the har- 
monic analysis of the corresponding Hilbcrt space C 2 (ji), i.e., the usual completion, 
now defined from 



(i.2i) ml- / i/i 2 *. 
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The simplest examples of this would be the following Cantor constructions: 
d=l, R = 4, = I = [0, 1] the unit-interval, and B = {0, 2}; 



or, alternatively, 



R = 3 and £ = {0,2}. 



The corresponding fractal dimensions are 1/2 = In2/ln4, respectively, 0.631 . . . = 
In 2/ In 3. As we show in [JoPe5-7], there are sets A C 1 such that for the first 
example, C? (/i) has a harmonic analysis which arises as a limit of the spectral pairs 
indicated in the construction (1.15)— (1.16) above, but this construction fails for 
the second example (the ternary Cantor set), and we will see in this paper that 
there is a limit construction which is still analogous, but it is based instead on the 
generalized spectral pairs, see (1.5) above. 

Following (1.3), we may define a transform jj, for Borel probability measures \x 
on R d : 

fi(X) := J e\(x) dn{x) 

and 

(1.22) := J exfdfi 

as functions in M. d . In his recent papers [Strl-2], Strichartz considers integrals on 
the form 



(1.23) I fd^(s) 



ds 



where a is a number (typically a fraction related to the fractional dimension of /z) , 
Ej a family of sets, Ej C Ej + i, such that < m d (Ej) < oo, and [J . Ej =R d . He 

shows that a may be chosen such that if C and C denote the respective lim inf j and 
limsupj in (1.23), then there are positive constants k, K such that 

(1-24) k\\f\\l <£<£< K\\f\\l. 

Compare this to the frame property (1.5) above. We shall be interested here in 
discrete versions of these spectral estimates; specifically, we shall give conditions 
on the dual affine system, as in (1.8)-(1.9), which guarantees an estimate as that 
of Strichartz in (1.23), but now with the lim inf and lim sup's defined instead from 
the sequence 

(1-25) (jf[y ]T |(7^)(A) 

The transform under the summation is defined in (1.22) and is the same as the one 
used by Strichartz. But neither there, nor here, is it clear when the j —> oo limit 
exists. In our setting, we illustrate this issue with results and examples in Section 
3 below. 
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Frame estimates as in (1.5) play a role in the theory of wavelets, and fractal 
geometry which is analogous to that played by the classical Garding inequalities 
in the theory of elliptic partial differential equations (p.d.e.). A main result of 
ours (details in Sections 2 and 3 below) states that, in the present setting, there 
is a finite-dimensional matrix A = M*M such that a given positive lower spectral 
bound for A implies an a priori "smoothness" for our fractal measure /z, and the 
corresponding sequence of transforms (1.22) which is the counterpart of a classical 
p.d.e. -Garding estimate. Specifically we show that a lower bound on the spectrum 
of A implies the following a priori estimate (for some e G M + ) 

(1.26) liminf (l||y ]T \(fd£) (\)\* > e J |/| 2 rf M 

for our fractal measure n (= lim^oo fij in the Hutchinson metric), which in turn 
arises as a solution to 

(1-27) M=^|E^°^ 1 

1 1 beB 

within the probability measures on R d . The equation (1.27) plays a role here which 
is quite analogous to that of an elliptic p.d.e. in the classical case. This is perhaps 
more clear if (1.27) is viewed as a special case of a more general class of difference 
equations, viz., 

(1.28) jj,= ^w b (x) fio o-^ 1 

beB 

where {u>b{-)}beB is a given family of functions on M. d . For more on this, see also 
[Str2]. The interpretation of (1.28) is the identity 

I fdn = Y, I w b (a b (x))f(a b (x))dfi(x) (for V/ G C c (R d )). 
J beB J 

Our dual affine systems are said to be elliptic if the two finite translation sets 
B and L determine a matrix A = M*M with (positive) spectrum bounded away 
from zero. Our analogy to the Garding case is the result (with suitable technical 
conditions) that ellipticity implies the a priori estimate (1.26). 

2. The Initial Step for the Iteration 

2.1 Pairs in R d . Let d be fixed, d — 1,2,..., and let m,d be Lebesgue measure. 
In Section 1 we introduced pairs of subsets (ft, A) of R d such that ft C K d is 
measurable and < m^Sl) < oo, and the set A serves as a spectrum for £1 in 
one of the two senses that we described; see especially (1.5) which corresponds to 
(fi, A) being a generalized spectral pair. If d = 1, then, in applications, f2 may be 
time and A frequency. From quantum mechanics, we may have d > 1, and points 
x = [x\ , • • • , Xd) in SI describe the position of a particle or system with many degrees 
of freedom, whereas points A = (Ai, . . . , Ad) in A describe the corresponding dual 
vector of momentum variables. 
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In (1.10) above, we considered the step from one set CI to a second set, CI + B, 
which is the union of translates of CI itself. The simplest special case of this would 
be when B = {0, a}, and the new set in x-space is, CI U (Cl + a), assuming again 
that 



Here Cl + a = {x + a : x <E Cl}. 

We shall be interested in generating new generalized spectral pairs from old ones, 
and eventually in the iteration of the process, leading thereby to affine fractals. The 
aim is to study the interplay between the original configuration in M. d on the one 
hand, and the iteration limit on the other hand. The original configuration has 
a harmonic analysis based on R d as an additive group and Lebesgue measure as 
Haar measure. The iteration fractal will have the canonical measure /i which is our 
substitute for Haar measure. 

One reason for the generalized spectral pairs is that, even if the starting point 
is a spectral pair, then the first, or one of the later affine iteration steps, typically 
takes us outside the restricted class of spectral pairs, but still we stay within the 
generalized class (i.e., g.s.p.). The generalized spectral pairs (Cl, A) are determined 
by the corresponding frame operator, F — F(n,\) from (1.3) and (1.7), i.e., the 
operator, F : f i— > (f(X))\ e \ of C 2 (Cl) into £ 2 (A); and we will now study how this 
operator F changes when we pass from SI to 



where a e R d is fixed, but chosen subject to (2.1). Recall the operator F maps 
from £ 2 (S1) into £ 2 (A). It is bounded, and it has zero kernel, i.e., the only solution 
to, Ff = 0, / G £ 2 (S1), is / = 0. It follows then from Hilbert space theory [BeFr] 
that F*F is a self-adjoint operator in £ 2 (S1). It has a spectral resolution 



where E(-) is an orthogonal projection-valued measure on R+. We will study the 
above mentioned problem by checking how (2.3) changes when we pass from SI 
to fli, where again Cli is determined as in (2.2) above. But we have in mind 
also more general steps, Cl i— > Cli, and we will follow through and supplying an 
accompanying modification, A ^ Ai, i.e., we will identify Ai such that (fii,Ai) 
is again a generalized spectral pair with an associated selfadjoint frame operator 
F*Fi, now in £ 2 (Sli) relative to Lebesgue measure on fli. The simplest special 
case of this construction is (see [JoPe7]) the case d = 1, O = [0, 1] = J, a € R, 
a > 1, and A = Z + /3, for some fixed (3 C R. It is clear that then (CI, A) is a spectral 
pair. This is just the classical Fourier transform. 

Example 2.1.1. We considered in [JoPe7] a > 1 for the non-overlap property 
(2.1), and we defined Cli and Ai as follows: 



(2.1) 



Cir\(Cl + a) = 0. 



(2.2) 



fti := f!U (CI + a) 



(2.3) 




(2.4) 



Oi := I U (I + a) = [0, 1] U [a, a + 1] 



and 



Ai := (Z + /3)Uo _1 Z 
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and we showed that (f2i,Ai) is a generalized spectral pair iff a £ Q and ^ 
Z + a _1 Z. The operator F*F\ is still well defined also if a is irrational, but we 
showed that the lower bound in the spectrum is then 0, and therefore that the 
lower bound in the frame estimate (1.5) is then violated for £li = [0, 1] U [a, a + 1], 
when a is irrational. These facts will follow from our Theorem 2.1.2 below, which is 
the corresponding general R d case. We will also give a completely explicit formula 
for the spectral resolution of the modified operator F^F\. 

If (f2, A) is a given generalized spectral pair (g.s.p. for short) in M. d , then it is 
clear from (1.5) that, for every vector I £ M. d , the pair (fi, A + I) is also a g.s.p.; 
and, similarly, we may make a translation of f2 and not affect the g.s.p. property. 
It follows that the condition, £ A, may be added, and the g.s.p. property is 
unaffected. In the case £ A, we shall need the dual discrete set A° from [JoPe3] 
and Lemma 1.1 above. We shall also need the similar duality (2.5) below, but now 
defined over the rationals, rather than the integers, i.e., 

(2.5) A°(Q) := {s £ M. d : s ■ A £ Q, VA £ A}. 

Theorem 2.1.2. Let (fi, A) and (f2 2 ,A 2 ) be two given spectral pairs in R d and 
assume that is in both A and A 2 and C fi 2 . Consider points a and in 
R d such that a satisfies the disjointness property (2.1), and now form a third pair 
(fli, Ai) as follows: 

(2.6) fii:=fiU(fi + a) and A 1 := (A + 0) U A 2 . 

If a £ A 2 , then the new pair (fii, Ai) is a g.s.p. if and only if 

(i) a i (A + 0)° and 

(ii) a£ A°(Q). 

In any case the spectrum of F±F\ is the set {r±(X) : A £ A} where 

(2.7) r±(A) = 2 + a ± yj(2 + a) 2 - 4a (1 - cos2tt(/3 + A) • a) 
and 

(2.8) a = 2m d (n 2 )/m d (ri). 

Proof. Let r be the period-2 automorphism of f2i which sends x £ il to x + a, and 
x £ Q + a to x — a; and define 

(2.9) Tf := for for all / £ £ 2 (^i). 
This operator is unitary in £ 2 (Qi) with spectral subspaces 

(2.10) H± := {/ £ £ 2 (fli) : Tf = ±/} 
and 



(2.11) 
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For elements in / in £ 2 (Sli) we may consider the transform / defined in (1.3), and 
the integration being over Oi, but we shall also need the restriction of / to O C ill 
and the corresponding transform, now with integration only over SI. The latter will 
be denoted /(•), specifically 

(2.12) /(A) := f ~e\(xj f(x)dx, for A G R d . 

Jn 

We shall be using the following shorthand notation 

(2.13) (X,x) := e x (x) = e airX ' x 

defined for all A, x G R d . When f± G H± (respectively) then the f± transforms 
compute with integration over fii = SI U (Si + a) as follows 

(2.14) /+(A) = (1 + <A^>)/+(A) 
and 

(2.15) /_(A) = (1-<AX)/_(A) 
for all A e M. d . Since a £ we conclude that 

/-(A 2 ) = forallA 2 eA 2 . 

From (i) in the theorem, we conclude that the two terms in Ai = (A + /3) U A 2 are 
disjoint, so the Ai-summation in 

(2-16) J2 |/( A i)f 

AiGAi 

may be calculated separately for A G A + f3, and then for A G A 2 . If the two sets 
in the union Ai := (A + [3) U A 2 arc not disjoint, we will still get a frame estimate, 
but then in the following modified form, 

\s< £ |/(A)| 2 <s 

AGAi 

where s := Ea £ a+/3 l/( A )| 2 + EagA 2 I/( a )| 2 5 but thc disjointness is implied by (i) 
and a G A^. With the H±-decomposition, we get 

(2.17) /(A) = (1 + (AT^))/ + (A) + (1-(AT^))/_(A) for all A G A + 0; 
and 



(2.18) /(A 2 ) = 2/+(A 2 ) for all A 2 G A 2 . 
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It follows that 



E I/m =4 e \m^) 



a 2 gA 2 



AigA 2 



4m d (fi 2 ) x / |/+| | 2 dx 
Jn 2 



= 4m d (tt 2 ) x 
, m d (Q, 2 ) 

m d (n) 



= 4- 



|/+| 2 ^ 



E 



AGA+/3 



where we used that 51 C f^2 and the assumption that both (f2, A) and ($72, A2) are 
given to be spectral pairs. But if they are known instead to be only g.s.p.'s, then 
the same calculation would give the estimate 



(2.19) 



E |/(*2) >p E |/+( A ) 



A 2 eA 2 



AGA+/3 



for some positive constant p (depending on the respective frame constants for the 
two g.s.p.'s). It will be convenient to write it in the form p = 2a, but this a will 
not be the one from (2.8) in the theorem. 

Returning to the summation (2.16), and still recalling the decomposition / = 
/+ + /_ from (2.11), we get 



E |/(Ai) =2a E M 

AiGAi AGA+/3 



E \(l + (\,a))x x + {l-(\,a))y> 

AGA+/3 



where the following shorthand notation was introduced: 

(2.20) x x ■= /+(A) and y x := /_(A) for A G A + f3. 

But this is the quadratic form which for each A e A + (3 is given by the 2 by 2 
matrix 



(2.21) 



2 ■ 



a + 1 + cos(27rA • a) i sin(27rA • a) 
—i sin(2-7rA - a) 1 — cos(27rA • a) 



and which by simple algebra is unitarily equivalent to the diagonal one 
where the roots r± are given by 



r+ 




r± = 2 + a ± x/(2 + a) 2 - 4a(l - cos(2ttA • a)) 

and this is formula (2.7) where points in A + /3 are written in the form A + /3 (A G A). 

If (i)-(h) both hold, it is clear that r±(A) > for all A e A, see formula (2.7) in 
the theorem, and (on account of (ii)) as A varies in A, the set {r±(A) : A e A} will 
then in fact be finite. Writing out the quadratic form (2.16)-(2.17), and using 

(/,F 1 *F 1 /)= E l/(Ai)| 2 
AeAi 
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for the sum in (2.16) above, we conclude that there is an orthogonal decomposition 
of £ 2 (fii) with corresponding projections {-E±(A)}agA such that 

(2.22) F*F, = J2 r+(X)E+(X) + £ r_(A)£L(A). 

AGA AgA 

Notice that each of the non-zero projections (in the rational case) 

J2{E±(X) : X e A, r ± (A)=£±} 

A 

is infinite-dimensional. Moreover, the rationality condition (ii) guarantees that the 
numbers {r_(A) : A e A} do not accumulate at 0. This means that we get the 
lower estimate in condition (1.5) for (fii,Ai) being a generalized spectral pair. 
Conversely, if is in the closure of {r_(A) : A e A} for some (fixed) a and (3, then 
this means that the lower bound estimate (1.5) for (fii, Ai) fails; and, by (2.7) and 
(2.21), this means that one of the conditions (i), or (ii), must be violated. □ 

We also get from the proof the following: 

Corollary 2.1.3. Let the setup be as in the theorem, but assume instead that (fi, A) 
and (fi 2 , A2) are only g.s.p. 's. Then we still get that the g.s.p. property for (fii, Ai) 
from (2.1) and (2.6) is equivalent to conditions (i)-(ii), but we will then not have 
the exact formula for the frame operator F*Fi in (2.3), but instead of (2.22) only 
the estimate 

(2.23) F*Fi > ^(r + (A) J B+(A) + r_(A)S_(A)) 

AGA 

and the number a in (2.19) and (2.7) will depend on the given frame estimates for 
the two g.s.p. 's (0,A) and (fi 2 ,A 2 ). 

Proof. We have already read off the lower bounds from the proof above. The upper 
bound for the new pair (fii, Ai) in (2.2) and (2.6) may be gotten from an application 
of [JoPc7; Theorem 4.2]. □ 

In summary, the difference between the respective cases, spectral pair and g.s.p., 
amounts to the difference between the equality in (2.22) and the estimate in (2.23). 

Remark 2.1.4- In Theorem 2.1.2 we construct new generalized spectral pairs 
(fii,Ai) from two given spectral pairs (fi,A) and (fi2,A 2 ). We want to choose 
A and a finite subset B C R d such that OsB, the non-overlapping property (1.12) 
(generalizing (2.1)) holds for (fi, .B), and B C A£. Moreover fi 2 is extending fi, 
i.e., we assume the inclusion fi C fi 2 . The theorem is about B = {0, a} where a is 
a single vector in M. d , and here we point out the modifications which are entailed 
by the generalization to when B contains more than two points, i.e., \B\ > 2. The 
above proof shows that when the frame operator 

(2.24) Fi : C 2 (Q + B) ( 2 {Ki) 

(where Ai = A U A 2 and A is some spectrum for fi) is introduced, we then still 
have the formula (2.22) above for F*F\ and the projections E±(X) are indexed by 
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A, and they are projections in £ 2 B . However the eigenvalues {r±(A) : A 6 Aj will 
be given by a different formula (see details later). But it follows that we cannot 
have the spectrum of F*F\ bounded away from 0. Specifically (2.22) amounts to 
the decomposition 

(2.25) F * Fl = 0Q (A ) 

AeA 

where each Q(X) : £ 2 B — > £ 2 B is a selfadjoint operator of rank 2. Since &\m.£ B = 
\B\ > 2, it follows that Q(A) cannot be strictly positive, and therefore, by (2.21), 
that must be in the spectrum of F*F\. This means, by (2.7), that we cannot 
have a positive lower bound in the frame estimate (1.5) for F\. 

The argument for why each Q(X) must be of rank < 2 for each A e A is based on 
a modification of the quadratic form (2.25) in the proof of Theorem 2.1.2. When 
B is given, introduce the function on R d : 

(2-26) <M0~£M. 

bEB 

and note that Q(X) may be represented as a positive definite function on ~R d x M. d , 
or rather on B* x B* where B* is a subset of R d with \B*\ = \B\, the positive 
definite function being 

(2.27) Q\(£, O = ^b(I)vb(O^St + <Pb{Z-\)<Pb{Il' A) 
for all (£,£') EB* xB*. 

2.2 Translation and scaling. In Section 1 we outlined the construction of new 
generalized spectral pairs (g.s.p.) on the form (£l + B, A + L) where (f2, A) is a given 
g.s.p. in R d , and B and L arc finite subsets of IR'' such that we have non-overlap 
of distinct translates f2 + b as b varies over B, and similarly for {A + £ : £ e L}. In 
Section 2.1 we considered a different candidate for a possible spectrum of tt + B, but 
in this section we will give conditions for getting a spectrum of the form A + L := 
{J £eL (A + £), again with non-overlap for the distinct ^-translates. 

Theorem 2.2.1. Let (fi, A) be a generalized spectral pair (g.s.p.) in R d with frame 
constants k, K as defined in (1.5), and let B, L be two finite subsets in M. d satisfying 
the respective non-overlapping conditions and suppose 

(2.28) B c A°. 

Let M. := ((b,£)) be the corresponding matrix where (b,£) := e aixh ' i are the en- 
tries, rows indexed by L and columns by B. Suppose A'1*A / 1 is bounded below, or 
equivalently that it is invertible in £ 2 B . With constants p, P such that 

(2.29) p\\v\\p B < (v,M*Mv) el < P\\v\\ t% forVvE£ 2 B , 

we have p > 0, and (fl+B, A+L) is again a g.s.p. with frame constants multiplying, 
i.e., the respective lower and upper frame constants, k\ respectively, K\ are given 
by 



(2.30) 



fci = kp and K\ = KP. 
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Proof. Consider four constants A\, A 2 , k, and K such that < k < K. We will be 
interested in double estimates on the form 

(2.31) kA 2 <Ai< KA 2 , 
and we will summarize this in the notation 

(2.32) A 1 = (k,K)A 2 . 

From the assumptions on the g.s.p. (Q, A) we get (2.32) satisfied for 

A6A 

and 



A 2 = I \f(x)\ 2 dx=\\ff c2(ny 



I l/( 
Jn 

The object is to establish that if A[ := J2a+l l/( A )| 2 > 



2 



A' 2 := / \f{x)\ 2 dx = ||., , k ., ( , H Li ... 
Jn+B 

and / eC 2 (Q + B), then 

(2.33) A'^ikuK^A'z 

where the new constants are given in (2.30). When the Fourier transform is on 
O + B it is denoted /, and when it is calculated on for the restriction, it is 
denoted /. Introducing the operator 

(2.34) T t f{x) = Y^f{x + b) 

bEB 

from C 2 {Q + B) to C 2 {Q) we then get 

f(\ + £)=fTf(X + £), 

and 



(2.35) Yl \ T 'f ( A + £ ) 2 - ^ X ) x Wf&w 

AeA 

for all / e C 2 (fl + B) and all I G L. Recall that each I + A serves as a spectrum 
for f2 and the frame bound is the same. 

For the second summation we get, using assumption (2.29), 

E / iwi 2 ^ = EEE(^- fc '> / W+b)f(x+b')dx 

e Jn b b 1 e 

= (p,p) [ Ei/^+^i 2 ^ 

Jn b 

= (P,P)[ l/l 2 
Jn+B 



12 dx 



in+B 

1 1 2 



(p,P)\\f\\Un + B), 
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and the result now follows by substitution. □ 

Remark. Systems of operators like (2.34) abound in wave-packet analysis, and in 
wavelet theory, see e.g., [JoPe3-6] and [Mey]. The corresponding adjoint operators, 
T; : C 2 (Sl) -> C 2 (Sl + B), are given by 

T*tf(x) = Xn+b(x)(£, b)f(x - b) for V/ G C 2 (Sl), 

b£B 

and it follows that 

7>T; = if B (£' - £)I C 2 {n) for W, £' G L, 
where <pb is given in (2.26), i.e., 

mo ■■= E< fo ^> = E ei27rK for v ^ e Rd - 

beB bEB 

In the special case where ipb(£' — £) = \B\5gj', we then get an associated orthogonal 
decomposition of C 2 (Sl+B), corresponding to the disjoint decomposition, [J beB (Q+ 
b) , into translates of a fixed domain ft (null-overlap of distinct translates is allowed) . 
A main point in the present paper is to show that wave-packet analysis carries over 
to the almost orthogonal case. 

The next result is about stability of the exactness property (see Section 1) in the 
step 

(2.36) (SI, A) (SI + B, A + L) 

in the category of g.s.p. Recall exactness for a given g.s.p. (SI, A) is the property 
that the associated operator F : C 2 (Sl) — > £ 2 (A) maps onto £ 2 (A). So it is equivalent 
to the adjoint 

F* : £ 2 (K) -» C 2 (Sl) 
having zero kernel, see (1.7) above. 

Corollary 2.2.2. If the initial g.s.p. (Sl,K) in the construction (2.36) from The- 
orem 2.2.1 is exact, and if in addition to B C A°, it is assumed that the matrix 
M. = ((b,£)) is invertible, i.e., \B\ = \L\ and M. invertible as a square matrix, then 
it follows that the new g.s.p. (SI + B, A + L) is also exact. 

Proof (sketch). Let (£) e £ 2 (A + L) and assume that F*£ = in C 2 (Sl + B) where 
Fi is the frame operator for the new g.s.p. in the Theorem 2.2.1 construction. Then 
the double sum 

eeL xeA 

must vanish a.e. x G SI, and for all b € B. Since B C A°, we get 

l£L AGA 

Invertibility of M yields 

]r& +A <^+A,x) = o 



A SPECTRUM FOR AFFINE FRACTALS 



15 



for all I € L and a.e. iefl. Using exactness of (Cl, A) we conclude that £g + \ = 
for all £ G L and all A 6 A. □ 

The final result is about the multiplicative property. Let R be a d by d matrix 
over the reals and assume that all its eigenvalues satisfy |Aj| > 1. The transpose 
is denoted R*. For the lemma, we just need R to be invertiblc, but the stronger 
property will be needed in the applications in the next section. 

Lemma 2.2.3. Suppose (Cl, A) is a g.s.p. in M. d and R is as specified. Then 
(R~ 1 Cl,R*A) is also a g.s.p. and the frame constants (k,K) are the same for the 
two g.s.p. 's. 

Proof. This simple argument is based on the transform formula for the Lebesgue 
measure in M. d and we leave details to the reader. □ 

2.3 Reversing the construction. We have seen how to get new generalized 
spectral pairs from old ones by the step Cl ^ Cl + B where Cl has a spectrum in the 
generalized sense and when B is a finite set such that the translates {Cl + b : b G B} 
do not overlap. In this section we show that the process can be reversed but this 
reversion is surprisingly subtle, and it is relevant in tiling theory, see [JoPe3-4] and 
especially [LaWa] . We show that if B has the non-overlap property and if there is 
a set A such that (Cl + B, A) is a generalized spectral pair, then there is an explicit 
construction of a second set Aq such that (Cl, An) is then also a generalized spectral 
pair. But as pointed out in [LaWa] even if d = 1 and Cl + B is a Z fundamental 
domain, i.e., if (Cl + B, Z) is a spectral pair (in the strict sense) in R, then typically 
(fi, Aq) will still only be a generalized spectral pair. 

Theorem 2.3.1. 

(i) Let Cl, A, B be subsets of R d , with B a finite set, and G A. Suppose 
B+fl = {Ji, e g(b+£T) is non-overlapping, and that (B+Ct, A) is a generalized 
spectral pair. Assume (e , e\)n+B = for all X e A\{0}. Let 



Then (O, An) is a generalized spectral pair. 
(ii) If (B + Q,A) has frame constants (k,K), then (O, An) has frame constants 



(k/\B\,K). 

Proof. For functions g on R d , let Z(g) = {A e A : g(X) = 0}, then A f2 = Z(xn)U(0). 
Let A e R d , then 



Hence, A\(0) = Z(x B+n ) = Z(ip B ) U Z(xa); in particular <^ S (A) = for all 
A G A\A n . Let / G C 2 (fl), and define g G C 2 (B + CI) by g(b + x) = f(x) for b G B, 



A Q :={AGA: X n(A)=0}u{0}. 
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x e 0, (i.e., g — 



To/ in the notation used above). 



k\B\\\f\\n 



k\\9\\ 2 B+ n<J2\9W\ 2 



agA 



2 



E J] / (&,A)(x,A)/(x)dx 



AgA bG-B 



E <MA)/(A) 



E Vb(A)/(A) 



AGAfi 



< 



\b\ 2 E i/»i 2 - 



AGAs-2 



This gives the lower bound. The upper bound follows from An C A. □ 

There is a much stronger version of the assumptions in Theorem 2.3.1 which 
stills gives an interesting application to the classification of spectral pairs in R d . 
First note that a subset V C R d is a fundamental domain if it is measurable and 
there is a lattice T in M. d such that D + T = R d with non-overlapping translates 
{V + 7 : 7 G T}. We say that (£>, T) is a lattice tiling. In [JoPe3] we considered 
subsets in R d on the form T> + B with B finite and non-overlap for the translates 
{V + b : b G B}. A main question is to decide the structure of spectral pairs on 
the form (D + B, A) when the sets V and B are specified as stated. (Note that if 
B cT, then the non-overlap property will be automatic.) If then G A we clearly 
have the assumption, 



satisfied. But if, as in Theorem 2.3.1, {T> + B,h) is only assumed to be a g.s.p., 
then this may not be the case. In Example 2.1.1 above, (2.37) is satisfied in fact iff 
a = 1 which is really the degenerate case. The following result serves as a converse 
to [JoPe3; Theorem 6.1] and classifies an important subclass of spectral pairs in M. d . 

Theorem 2.3.2. Let T> be a fundamental domain in M. d , and let A, B be subsets in 
M. d such that B is finite and satisfies the non-overlap property relative to V, G A, 
and (V + B,A) is a spectral pair in M. d . Suppose that there is some lattice T such 
that (V, T) is a lattice tiling, and that we have the following inclusion 



Then it follows that there is a second finite subset L C K d of same cardinality as B 
such that A = L + T° where T° denotes the lattice which is dual to T. 

Proof. From Theorem 2.3.1 we conclude that (D, Ax>) is a g.s.p. where 



(2.37) 





{A G A : xv(\) = 0} G r°. 



(2.38) 



Ax, = {AgA:xd(A) = 0}U{0}. 



But since (V, T) is a lattice tiling, and the pair (V + B, A) is a spectral pair, we 
conclude for the set A v in (2.38) that A v C T°nA. So A v C Y° and {e s \ v : s G T°} 
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is an orthogonal basis for C 2 (T>). Since {e s |x> : s G Ax>} is total in C 2 (T>), (recall 
that the pair (D, Ax>) was a g.s.p. by Theorem 2.3.1), we conclude that 



and so T° C A. If £\ G A\T° we may apply the same argument and Theorem 2.3.1 
to the spectral pair (D + B, A — £\) and conclude that r° C A — £\, or equivalcntly 
{0,4} + T° C A. If £ 2 G A\({0,£i} + then the same argument would give 
{0, £i, £2} + ^° C A. But an application of [JoPe3; Theorem 6.1] and [Lan; Theorem 
4] shows that the inductive process must stop, and then the corresponding matrix 
(\ i B\~ 1 / 2 e' l27 ' l] ' bk ) . k will be unitary. So in particular the set L :— {0,£j} will have 
the same cardinality as B, and L + T° = A. □ 

Further details and applications of this result to the classification problem for 
spectral pairs in R d will be given in a subsequent paper by the coauthors. 

2.4 Spectral pairs in one dimension. This section serves as an appendix to 
Section 2.3 above, and shows that, when d = 1, and when the fundamental domain 
T> from Theorem 2.3.2 is further assumed to be the unit-interval / = (0,1), then 
we have all the extra assumptions in Theorem 2.3.2 automatically satisfied. As a 
corollary, we therefore get a classification of all spectral pairs (for d = 1) on the 
form (f2, A) where f2 is a finite (non-overlapping) union of translates of /. Specifics 
are worked out in detail for the convenience of the reader. 

The next result, as the one before it, has a converse, which in fact follows from 
our earlier paper [JoPe3; Sections 6-8]. Putting the results from the two papers 
together, we then get a complete classification of all spectral pairs (fi,A) in one 
dimension where the set O is built from translates of an interval; specifically, Q = 
I + B = \J beB (I + b), with non-overlap, where / = (0, 1), and B C M is a finite 
subset. 

Corollary 2.4.1. Let I = (0, 1), B C R be finite, B + I be non- overlapping, and 
let ACM. Assume that (B + I, A) is a spectral pair. Then A = L + Z for some 
finite set L, \L\ = \B\, L + Z = {J f eL £ + Z is non-overlapping, and b — b' G Z for 
all b, b' C B. Moreover, the matrix (|_B| _1 / 2 (6, £))beB,£eL must be unitary. 

Proof. Let I e A, and define 



Hence <a{t) = iff t e Z\(0), and A/ C Z. By Theorem 2.3.1, (J, A/) is a 
generalized spectral pair; in particular, {eA : A G A/} is total in C 2 (I). But 
{e\ : A G Z} is an ONB for £ 2 {I), and A/ C Z; thus A/ = Z. It follows that 
£ + Z = ^ + A/CA. Inductively, suppose {£\, . . . , £ p } + Z C A, and let £ G 
A\({£ 1, . . . , £ p } + Z). Arguing as above, we get {£, £\, . . . , £ p } + Z C A. By [Lan, 
Theorem 4] we have 



A p = r°nA = r° 



Aj :={\EA-£:xiW 



0}U{0}, 



where 




d + (A) = limsup ( r 1 max |A n B a (r)\ I < mi(O) 
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where B a (r) = [a,a + r]. Hence the construction above must terminate, and A = 
L + Z for some finite set LcA. By construction, £ — £' Z for all £, £' E L with 
£ ^ £'. Let 

pb(*) = £(M>, tem. 

bEB 

Then 

<p B (£-£' + z-z') -xi{£-£' + z-z') = f el^7e e+z = 

JB+I 

if £' + z' ^ £ + z. It follows that 

ip B (£-£' + z- z') = 0, for all z, z' &Z,£,£'eL 
with £ ^ £'. For z E Z, let M z be the |B| x \L\ matrix 

M z = ((b,£ + z)) beB ,eeL 

(bEB labelling the rows, and £ E L labelling the columns). Then M z has orthog- 
onal columns; hence \B\ > \L\. By totality of {ei +z : £ E L, z E Z} in C 2 (B + I), 
we have \L\ < \B\. Thus \B\ — \L\, and M z is a Hadamard matrix for all z E Z. 

For the rest of the proof, we will without loss of generality assume that E B 
and 6 I. Define vectors v n , w e E £ 2 B by v n = ((b,n)) beB , w t = {(b,£))beB for 
n E Z and £ E L. Then 

v n -w t = y^(b,£- n) = ip B (£ - n). 

b 

Hence v n ■ W£ — for £ E L\(Q), n € Z. The vectors in {w£ : £ E L\(0)} are 
mutually orthogonal, and hence span an (|L| — l)-dimensional subspace in £ B . It 
follows that there are scalars c n E C such that 

v n = c n v = c n w 

for all n E Z. (Note v = (1, 1, . . . , 1).) 

So (b,n) is independent of b E B. Therefore (b, 1} = (0, 1) = 1, and it follows 
that bEB; i.e., B C Z. □ 

3. The Spectrum of the Limiting Fractal 

3.1 Iteration. In this section we consider results and examples which realize limit 
properties for the approximating discrete transforms 

Fnf = I d^ n = J e\f d/j, n 

which enter into our iteration theory as outlined in Section 1 above. For a given 
generalized spectral pair (f2,A) in R d we shall need Lemma 1.1 about the discrete 
set K := A°; and we shall also need the simple observation that every translate 
(ft + a, A + {3) is again a g.s.p. with the same frame constants as the original one 
(SI, A). 
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Let (f2, A) be a generalized spectral pair, let B and L be finite subsets of R d , let 
R be an invertible d x d matrix with real entries. Define affine maps of R d , a\, and 
t £ , for b G B and £ G L by 

(3.1) cr 6 x = iT^x + fr) and T £ s = i?*(s + £) 

for all x, s G M d . Inductively define sequences (fl n ), (A„), (n — 0, 1, . . . ) of sets by 

(3.2) ft = n, A = A, 

(3.3) n n+1 = cr(n n ) = (J cr 6 (fi n ), 

fees 

and 

(3.4) A„ +1 = r(A„) = [J n{A n ). 

leL 

The following five conditions will be standing assumptions for our duality 
iteration-algorithm: 

(i) the matrix M = ((£,b))e e L,beB, where (£,t) — e^(t) — e* 27 ^'*, has rank \B\, 

(ii) R{K) cK, 

(iii) B C K , 

(iv) L C R—^K ), and 

(v) the sum L + K° is direct; i.e., {£ + K°) n + F) = for W ^ £' in L. 

Here -FT = A° as before. Then (Q„, A„), n = 0, 1, . . . , is a sequence of generalized 
spectral pairs, as will be shown below. In addition to the results from Section 2 
above we shall also need the following four lemmas: 

Lemma 3.1.1. Let ft, A, and R be given as described above, and assume that 
conditions (i)-(v) hold. Let (Q„) and (A„) be the iterations. We then have: For 
each n G {0, 1, 2, . . . } 

(a) BcA°. 

(b) <7(,(f2„) n <Tfe/(O n ) is a null-set whenever b, b' G B and b =/= b' . 

(c) n(A n ) n n>(A n ) = whenever £, £' G L and £ ^ £'. 

Proof. 

(a) Recall A„ = {ft*™ A + YT k =i R * k£ k : A G A, 4 G L}. Using (ii), K = A°, 
(iv), and (iii), we get R*K° C K° , A C A°° = fl*L C and 
K° C B°. Hence, A„ C F C B°. Taking polars in turn yields, B C B°° C 
ftc A°. 

(c) From the proof of (a), we have A„ C AT ; hence the desired conclusion 
follows from the asumption (v) that (K° + £) n (F + £') = 0. 

(b) Claim: If {e p : p G ft} is total in C 2 {S), then 5 n (5 + 77) is a null-set for 
any 77 G ft°\(0). Let f(x) = 1 if x G S(t?) = S n (S + 77), and /(x) = if 
x G S\S(t}). Note that 5(7?) C S, 5(77) - 77 C S, and e p (x) = e p (x - 77) for 
all x G 5(77) and all p G ft. But / (x) = 1^0 = /(x — 77) when x G 5(77), so 
the totality of {e p : p G ft} implies that 77 = 0, or 5(77) is a null-set. This 
proves the claim. 
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Completion of the proof of (b). We will show, by induction on n, that (fi n , A„) 
is a generalized spectral pair, and that ab(fl n ) H ay^n) is a null-set. For n = 0, 
this follows from the assumption that (f2, A) is a generalized spectral pair, and the 
claim. 

If (fi n , A„) is a generalized spectral pair and the intersections <Tb(Q n ) n ov(fi„) 
are null-sets, then it follows from part (a), part (c), and Theorem 2.2.1 and Lemma 
2.2.3, that (f2„ + i, A„_|_i) is a generalized spectral pair, but then the claim, together 
with the observation A n+ i C K (from the proof of part (a)), finishes the proof. □ 

Let A denote a Borel set in M. d . Let m(A) denote the Lebesgue measure of the 
set A. Using Lemma 3.1.1(b), we then note that 



(3.5) 0„ = | J] R^bk : b k G B j + R~ n Sl 

is without overlap (redundancy), so we may define a probability measure [i n on 
R d by 

I det R\ n 

(3.6) ^(A) = J L_^ md{An n n) 

for Lebesgue measurable sets A. 

Lemma 3.1.2. The measures [i n defined above satisfy 



(3-7) Mn+i = j4E /J ' ,0(7 ' 

/or n = 0,1,.... 

Proof. Using Lemma 3.1.1(b) we get 



/ 



Ideti?!^ 1 f 

fdM= w^m LJ {x)dx 



dctmn+1 E / /(*)* 



|B|»+im(n)— Jok(nB) 
|dcti?|" 



|B|"+ 1 m(fi) 



E / f( (J bx)dx 

beB 



= TojE / f 0(J bdHn. □ 

1 1 beB J 



If all the eigenvalues Aj of i? have |Aj| > 1, then Hutchinson's theorem [Hut] tells 
us that the sequence (/z n ) converges to some limit probability measure fx satisfying 
t 1 = ]b~\ Sfc t lo<J b ' or equivalently, J f d\i — ^2 b J fo<j b d/j for all Borel functions 
/; and we now turn to the analysis of /i. 

For future reference we first record the simple 
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Lemma 3.1.3. If the constants a < A are such that 
(3-8) a\\ff n <£|/(A)| 2 <A||/|| 



for all f e C 2 (n), then 



I Mo 



/or ^/e£ 2 ( Mo ). 

From now on we will assume that a triple, A, B, L C R d , a sequence of probability 
measures /z n on M d , and an invertible <i x d matrix i? are given. We will suppose 
that the five conditions (i)-(v) are satisfied, and that 

(3-10) Hn+l = T^T ° °l ^ 

1 1 6GB 

for n = 0, 1, . . . . Notice that Lemma 3.1. 1(a), (c) remain valid, with the same proof. 

If all eigenvalues of R have modulus > 1, the theorem in [Hut] tells us that [i n 
converges to a probability measure fi satisfying (j, = jj^ J2beB M ° a b X - 

Lemma 3.1.4. For f e C 2 (fi n ), let 

;F„/(A) = J f(x){x, A) dfJk n (x). 
Then for f e £ 2 (^„ +1 ), I e L, A c A n , we have the transformation rule 

V b£B ) 

Recall that (x, A) = e x (x) = e i27lX - x . 
Proof. 

{F n+ if){T/\) = j f(x) (x,n\) d(j,n+i(x) 

= i4E / f( a b x ) (o- b x, TfX) dfi n (x) 



T5T$^(M) / (x,t)f(o- b x) (x,X)dfj, n (x) 



as desired. We used that (6, A) = 1; this follows from Lemma 3.1.1(a). □ 
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Theorem 3.1.5. Let A, L, B, R, fi n be as defined above. Assume there are 
constants k < K such that 



(3.11) 



fc||/||* < £ \T f(X)\ 2 <K\\ 



AeA 



for all f E £ 2 (^o)- Let M = ((b 7 £))i,£BJeL, an d let d < D be such that the \B\ by 
\B\ matrix M*M satisfies 



d\\v\\l m < (v,M*Mv) cm <D\\v\\ 2 clBl 



for all v e CH B I . Then 

d \ n 



(3.12) 



k 



\B\ 



< \Fnf(\)\ 2 <K 



AGA„ 



D_\ n 

\B\ 



for all f E C 2 {[i n ) and each n = 0, 1, 

Here, Ao = A and A„ + i = t(A„), as before, and T n is as in Lemma 3.1.4. Notice, 
by Lemma 3.1.3, that if (f2,A) is a generalized spectral pair and \xq 



(meaning /xo(A) = m ^ m(An^) for all Borel sets A). Then the hypothesis of this 
theorem is satisfied. 

Proof. The proof is by induction on n. Assume that 



fc (^)"ii/iiL^Ei^/wi 2 ^(|i)" 



for allfeC 2 (p n ). 

For / G £ 2 (tin+i), let 



then 



Tif = e t -^Y / (bJ)f°Vb, 



b£B 



EimL 



—y 

\B\ 2 f-^ 



b£B 



^EE(/°^E< fo ^>/°^') 

1 1 beBb'eB \ e I Mri 



1 



(v(x), M* Mv(x)) c \ B \ djjL n {x) 



where v(x) is the vector in C' 5 ' with entries (/ o ab(x))i, e B- Hence 
d 



\B\ 



• f \\v{x)\\l m d» n (x) <J2\\T e f\\ 



< 



W 2 



\v(x)\\ 2 :lBl dfj, n (x) 
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it follows that 

(*) ||rii/iiL +1 <Eii T ^iiL< i?tII/iiL +1 - 

11 ££L 11 



Now 



E l^n+l/(0| 2 =E E l*n+l/(T/A)| S 

?GA„ + i £eLAeA„ 

|2 



= EEi(«/)( A )i 5 

* A 

by Lemma 3.1.4, but then the inductive hypothesis yields 

k (wX J2\^f\\L ^ E i^ + i/(OI 2 <^(|t)"eii^/iiL- 

Combining this with (*) gives the desired frame estimate. □ 

3.2 Limiting estimates. Assume now expansitivity, i.e., that the eigenvalues of 
R all have modulus > 1. Then, as noted earlier, the sequence (fi n ) converges to 
a certain limit measure fi, so that ||/||^ n converges to ||/||^ for all continuous /. 
We could then ask for the existence of a renormalization constant T such that the 
terms 

(3.13) T n l^/(A)| 2 

AeA„ 

are comparable to ||/||^ n for all n. 

The next three results present a solution to this circle of asymptotic problems. 

Theorem 3.1.5 indicateds that this will be the case for T = provided d = D; 
i.e., provided 

(3.14) (v,M*Mv) ilBl =D\\v\\ 2 lBl 

I B I 

for all v 6 4 '. In this case the columns of M are orthogonal and D — \L\, the last 
fact used that all the entries in M have modulus 1. In particular, d = D if and 
only if £ £ei <& - V, £)=0 for all b, b' E B with b ji b'. 

Proposition 3.2.1. Let (CI, A) be a generalized spectral pair, with frame constants 
a < A. Suppose (i)-(v) are satisfied. If^ ieL {b — b',£) — for all b, b' e B with 
b^b', then 

(3.15) ^totII/IIm ^ liminf f Ifr)" E \r»fW\ 2 



an 



d 



(3.16) limsupfM") l^/(A)| 2 < 

for all f e C(E). 



^ „„,2 
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Corollary 3.2.2. If the original frame (CI, A) is tight (i.e., a = A), then 

A 



(3.17) i(f)"El^)! 2 - 



2 



aga„ md ^ 



In general fj, n — > fx as n — > oo, hence \\f\\^, n — > ||/||^, cwid -T 7 ™/ — * as n ^ oo 

for any continuous function f. 

Examples 3.2.3. We will use (fl,A) = ((0,1), Z) in the examples. Then 



I/II,C 2 (0,1) 

Aez 



for all / e C 2 (n): hence a = A = m(Q) = 1. Also K = K° = Z: hence the 
conditions (i)-(v) are: S C Z, i? G Z\(0), L C ^Z, if t, t 1 & L then i-i 1 and 
the matrix ((6, l))beB,eeL has rank We will always assume £ B, £ i. 

(i) Let B = {0, 2}, L = {0, ±}, R = 4. Then M = ((6, *)) = ^ J ^ ^ , hence 
d = D=\L\ = \B\ = 2, so that 



(3-18) ]T |^„/(A)| 



ll/ll* 

AeA„ 



as n — > oo, for any continuous /. Ao D Ai D A2 D • • • , and H^Lo 
C(L), where 



(3.19) C(L) = { E -ij^L, pG 

It is known (see [JoPe6]) that 

E l^/(A)| 2 =0 

Ae£(L) 

is possible for / ^ 0. 
(ii) Same as (i), except we replace L = {0, j} by L = {0, j, |, |}. Then A„ = A 
for all n. (The sequence /Lt„ is the same as the one in (i).) In this example 
Sagz l^>/WI 2 = +°° is possible, e.g., if f(x) = 1 all a;, then T^f = ft 
and it is not hard to show that ^ fi(2) = /t(2 • 4 J ) for any j e N. 

Now M= \ / I , so d = D = \L\ = 4, \B\ = 2, and therefore 




(3.20) Q)"El^™/( A )| 2 ^ 11/11' 

^ ' Aez 

as n — > cxo (for any continuous /). 
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(iii) Let B = {0, 2}, L = {0, 1/3}, R = 3. The M = ^ J gi i /3 ^ ; hence d = 1 



and D — 3, and therefore 



1 V / 3 

1 ||/IU < £ l(^/)(A)| 2 r ' ' 

AeA„ 



So in this case, renormalization appears to be impossible. In this example 
the limit measure /i is the ordinary Cantor measure on the middle-thirds 
Cantor-set. 

(iv) Same as (iii) except L = {0,1/3} is replaced by L = {0, |, |}. (So [i n 
and n are the ones from (iii), but A„ is changed.) Then A„ = A = Z for 

all n. Also M = [ 1 e li7T ? 3 ] : hence M*M = ( J J Y It follows that 

d = D = \L\ = 3; hence 



(3.21) (0"El^/( A )l 



2 

Aez 



as n — > exo, for all continuous /. 
(v) Suppose we pick L maximal, i.e., as a complete set of representatives for 
the cosets -j^Z/Z. If b e B, b ^ 0, then the rank condition prevents R 
from dividing b; that is, it means that a = e 4271 " 6 /^ ^ l. Hence X^tTc) 1 ° J ' = 
= 0, thus X*X = |L|7. We then get 



v 1 17 Aez 

for all continuous /. (See (ii) and (iv) above.) 



(3-22) Jim (tTt) £ l*»/(A)| 2 = ll/l 

agz 



2 



Acknowledgments 

The research going into the present paper was done while the coauthors partic- 
ipated in the two work conferences on wavelets and fractals in the Spring of 1994, 
one at the University of Pittsburgh (USA) in May, and the other at Finsterbergen 
(near Jena Universitat in Germany) in June, and we gratefully acknowledge gen- 
erous financial support from the respective sponsors of the two conferences. We 
also thank the individual organizers for stimulating meetings, especially Professors 
K.-S. Lau (University of Pittsburgh), and C. Bandt (Greifswald Universitat). We 
had helpful discussions with them, with the other participants, and especially with 
M. Lapidus, R. Strichartz, D. Mauldin, J. Lagarias, and Y. Wang. We are also 
grateful to the last mentioned two (Lagarias and Wang) for sending us in July 
a preprint version of their wonderful paper [LaWa] which inspired and motivated 
some of the present results. 

The work was supported in part by grants from the U.S. National Science Foun- 
dation, and from NATO. The first named author was also supported by a University 
of Iowa Faculty Scholarship. 



26 



PALLE E.T. JORGENSEN AND STEEN PEDERSEN 



References 

[Benl] J.J. Benedetto, Frame decompositions, sampling, and uncertainty principle inequalities, 

Wavelets (J.J. Benedetto et al., eds.), CRC Press, Boca Raton, 1994, pp. 247-304. 
[Ben2] , Stationary frames and spectral estimation, Probabilistic and stochastic methods 

in analysis (J.S. Byrnes, ed.), NATO 1991, Kluwcr, Dordrecht, 1992, pp. 117-161. 
[BeFr] J.J. Benedetto and M.W. Frazier, Introduction, Wavelets, Mathematics and Applications, 

CRC Press, Boca Raton, 1994. 
[BeTo] J.J. Benedetto and A. Teolis, Local frames, Mathematical Imaging: Wavelet Applications 

in Signal and Image Processing, SPFE, vol. 2034, 1993, pp. 310-321. 
[Bcu] A. Bcurling, Local harmonic analysis with some applications to differential operators, 

Some Recent Advances in the Basic Sciences, vol. 1, Academic Press, New York, 1966, 

pp. 109-125; The collected works of Arne Beurling, vol. 2, Birkhauser, Boston, 1989, 

pp. 299-315. 

[DuEa] R.J. Duffin and J.J. Eachus, Some notes on an expansion theorem of Paley and Wiener, 

Bull. Amcr. Math. Soc. 48 (1942), 850-855. 
[DuSc] R.J. Duffin and A. Schaeffer, A class of non-harmonic Fourier series, Trans. Amer. Math. 

Soc. 72 (1952), 341-366. 
[Fal] K.J. Falconer, The geometry of fractal sets, Cambridge University Press, Cambridge, 

1985. 

[Fug] B. Fuglcde, Commuting self-adjoint partial differential operators and a group theoretic 

problem, J. Funct. Anal. 16 (1974), 101-121. 
[Gab] J. -P. Gabardo, Weighted tight frames of exponentials on a finite interval, Monatsh. Math. 

116 (1993), 197-229. 

[Gro] K. Grochcnig, Describing functions: Atomic decompositions versus frames, Monatsh. 

Math. 112 (1991), 1-42. 
[Hut] J.E. Hutchinson, Fractals and self-similarity, Indiana Univ. Math. J. 30 (1981), 713-747. 
[Jol] P.E.T. Jorgensen, Spectral theory of finite-volume domains inW 1 , Adv. Math. 44 (1982), 

105-120. 

[Jo2] , Scattering theory for orthogonal wavelets, preprint, University of Iowa (1993). 

[Jo3] , Operators and representation theory: Canonical models for algebras of opera- 
tors arising in quantum mechanics, Notas de Mathematica, North-Holland, Amsterdam, 
December 1987. 

[JoPel] P.E.T. Jorgensen and S. Pedersen, Harmonic analysis and fractal limit-measures induced 
by representations of a certain C -algebra, J. Funct. Anal, (to appear). 

[JoPe2] , Harmonic analysis of fractal measures, Constr. Approx (to appear). 

[JoPe3] , Spectral theory for Borel sets in R n of finite measure, J. Funct. Anal. 107 

(1992), 72-104. 

[JoPe4] , Group theoretic and geometric properties of multivariable Fourier series, Expo- 
sition. Math. 12 (1993), no. 2, 1-24. 

[JoPe5] , Harmonic analysis of fractal measure induced by representations of a certain 

C* -algebra, Bull. Amcr. Math. Soc. 92 (1993), 228-234. 

[JoPe6] , Sur un probleme spectral algebrique, C.R. Acad. Sci. Paris Ser. I Math. 312 

(1991), 495-498. 

[JoPe7] , Local harmonic analysis for domains in R™ of finite measure, Preprint. 

[LaWa] J.C. Lagarias and Y. Wang, Tiling the line with one tile, Preprint. 

[Lan] H.J. Landau, Necessary density conditions for sampling and interpolation of certain 

entire functions, Acta Math. 117 (1967), 37-52. 
[Mey] Y. Meyer, Wavelets and operators, Different Perspectives on Wavelets (I. Daubcchics, 

ed.), Proc. Symp. Appl. Math., vol. 47, American Mathematical Society, Providence, RI, 

1993. 

[Olc] A.M. Olcvskii, Fourier series with respect to general orthogonal systems, Ergebn., vol. 

86, Springer- Verlag, New York, 1975. 
[PaWi] R.E.A.C. Paley and N. Wiener, Fourier transform in the complex domain, American 

Mathematical Society, Providence, RI, and New York, 1934. 
[Strl] R.S. Strichartz, Besicovitch meets Wiener: Fourier expansions and fractal measures, 

Bull. Amer. Math. Soc. 20 (1989), 55-59. 



A SPECTRUM FOR AFFINE FRACTALS 



[Str2] , Wavelets and selfaffine tilings, Constr. Approx. 9 (1993), 327-246. 

Department of Mathematics, University of Iowa, Iowa City, IA 52242, USA 
E-mail address: jorgenamath.uiowa.edu 

Department of Mathematics, Wright State University, Dayton, OH 45435 
E-mail address: spedersen@desire.wright.edu 



